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ABSTRACT 

Linear analysis of the formation of protostellar cores in planar magnetic interstellar clouds yields 
information about length scales involved in star formation. Combining these length scales with var- 
ious distributions of other environmental variables, (i.e., column density and mass-to-flux ratio) and 
applying Monte Carlo methods allow us to produce synthetic core mass functions (CMFs) for different 
environmental conditions. Our analysis shows that the shape of the CMF is directly dependent on the 
physical conditions of the cloud. Specifically, magnetic fields act to broaden the mass function and 
develop a high-mass tail while ambipolar diffusion will truncate this high-mass tail. In addition, we 
analyze the effect of small number statistics on the shape and high- mass slope of the synthetic CMFs. 
We find that observed core mass functions are severely statistically limited, which has a profound 
effect on the derived slope for the high-mass tail. 

Subject headings: diffusion - ISM: clouds - stars: formation - stars: luminosity function, mass function 
- ISM: magnetic fields - ISM: structure 



1. INTRODUCTION 

Observations of the stellar initial mass func- 
tion (IMF) and the core mass function (CMF) 
show similarities in the shape and high mass 
slope of these two functio ns (iMott^ et a l.l 1 19981; 
Testi & Sargent 1998; Johnstone et al. 2000; Alvcs ct al. 
2007): .Nutter k \y ard- Tho mpson. ,2007; .Simpson et al.. 
20081: lEnoch et al.l l2008: Sa davov et al.ll2010l among oth- 
ers). As such, much theoretical effort has been in- 
vested in order to explain these similarities. Vari- 
ous different approaches to this problem have been 
explored, including analytic and numerical studies 
which invo k e gravita t ional fragmentation or acc retion 
(fSiig [Togs'; 'Inutsukal !200ll; iBasu & Jones' '200i), tur- 
bulence (Padoan et al. 19 97; Padoan & Nordlund 2002; 
Ballesteros-Parcde s et al.l |2006; Hennebellc k. Chabrier 
2008 . 2009.) , independ ent stochastic processe s (.Larso: 



/UUo. jUUij), mdepena ent stocnastic processe s i.Larsoir 
1973 ; lElmegreenl 119971 ) and magnetic fields (|Dib et al.l 
2008( 1 ■ among others. Results of these studies vary from 
those which seem to agree with the fiducial Salpeter form, 
dN/dlogM (X M~" where a = 1.35 is the value of the 
Salpeter slope, to those that do not. 

The high mass s lope of the IMF was initially derived 
bv | Salpetei (Il955fl and later i mprov ed upon by iKroupal 
(|2002t ) and lChabrieri (|2008allhl . [2005h . Despite variations 
in observed and theoretically derived IMF slope values, 
it is often assumed that the shape and high mass slope 
of the IMF and CMF are identical and universal. From 
a theoretical view, such a one-to-one correspondence be- 
tween these two functions implies that high-mass cores 
beget high-mass stars and likewise for low-mass cores. 
The need for extensive simulations of how a complex of 
cores turns into a cluster of stars is simplified tremen- 
dously if it is assumed that each core will collapse into a 
single star with some mass loss to account for the mass 
shift between the CMF and IMF. 

The underlying tenet of universality is that all star- 
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forming regions are the same and undergo the same pro- 
cess to form stars, however recent observations and sim- 
ulations have started to reveal cracks in this assumption. 
In a study of the effect of turbulence on the form ation 
of the CMF, iHennebelle fc Chabrieii pOOSl . [20091) find 
that comparisons between their IMF and observations 
for different cloud conditions suggest that star formation 
should predominantly occur i n clou ds five times denser 
than characterized by ILarso"nl (jl98lD . This led them to 
question the universality of the IMF since, as they say, 
choosing different cloud parameters would lead to a dif- 
ferent CMF/IMF. Several recent studies of the IMF also 
tend to disagree with the assumed universality. Observa- 
tions of different star clusters in both the Milky Way and 
the Large Magellanic Clou d (LMC) show a wide scatter 
of slopes: a = 0.5 — 2.0 (jElmegreeiJ Il999[ ). A survey 
of high mass slope values for different stars (i.e., clus- 
ter stars versus association stars versus field stars) yields 
a wide range of values; a — 2.0 — 4.0 for extrerne fiel d 
stars to a = 1.5 — 2.0 for clu s ter st ars (jElmegreenll 1 99 71 ) . 
Further to this, lElmegreenl (|1999| ) shows that through 
stochastic fractal sampling of a cloud, the derived IMF 
slopes can vary fro m a as low as 1.0 to as high as 1.7. 
iClark et al.l (|2007D note that if the lifetime of a more 
massive core is longer than a less massive one, the slope 
of the CMF should be shallower i n order to obtain the 
IMF. Finally, iZaritskv et"all (j2012| ) show that there may 
be evidence for two distinct stellar IMFs that depend on 
the age and metallicity of the cluster in question. Based 
on the above evidence and arguments, it is not clear why 
one should insist on using a — 1.35 as the universal slope 
for both the CMF and IMF. 

The majority of the work in this area has focused on 
the effects of turbulence within the molecular clouds on 
the formation and shape of the CMF. Research which 
considers the effect of magnetic fields and ambipolar 
diffusion on the CMF is sparse. iKunz fc Mouschoviaj 
(■2009. ) used the results of a non-idea l MHD linear 
analysis of a partially ionized sheet (jMortonl 119911 : 
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iCiolek fc Basiill2006[ ) to generate a broad CMF, assum- 
ing ambipolar-diffusion initiated core formation. Their 
model assumed subcritical to critical initial conditions 
with a uniform distribution of mass-to-flux ratios be- 
tween 0.1 and 1.0 times the critical value for gravita- 
tional instability (see Section 5 for more discussion of 
their model). 

In this study, we use the results of the linear analysis of 
a partially ionized sheet along with a lognormal distribu- 
tion of initial column density and various distributions of 
mass-to-flux ratio. We explore both subcritical and su- 
percritical initial conditions. Mildly supercritical initial 
conditions are the most likely to l ead to massive core for - 
mation, as seen in e.g., Figure 2 of lCiolek fc Basul ()2006D . 
Furthermore, we use a lognormal distribution of column 
densities, as expected in molecular clouds on both the- 
oretic al grounds for a turbule nt medium f Padoan et al.' 
I1997I) and from observations (jKainulainen et al. 2009). 
The aim of this paper is two fold. In the first part 
we show the effects of a magnetic field on the shape of 
the CMF. Starting from an assumption of lognormal col- 
umn density probability we show the broadening effect 
of neutral-ion drift via ambipolar diffusion and differing 
mass-to-flux ratio distributions. In the second part, we 
address the inherent limitations of observed core mass 
functions, i.e. sample size and bin size. Specifically, we 
aim to compare small sample synthetic CMFs to large 
sample synthetic CMFs to show effect of small number 
statistics on the observed features of the CMF. In Section 
2 we outline our model and methods for constructing our 
synthetic CMFs. Section 3 shows the results for the dif- 
ferent distribution models considered. Section 4 shows 
the effect of small number statistics and the variance in 
derived analytic slopes. Finally Sections 5 and 6 give our 
discussion and conclusions. 

2. SYNTHETIC CORE MASS FUNCTIONS 

To better understand the effects of the environment 
on the shape and peak of the core mass function, we 
produce synthetic CMFs (synCMFs) based upon vary- 
ing physics and properties of molecular clouds. These 
include the column density ((T„,o), ionization fraction 
(Xi = log [?T-e /"■_?/]): mass-to-flux ratio (/^o): and neutral 
ion-collision time (T„i). The synCMFs are produced by 
randomly sampling predefined column density and mass- 
to-flux ratio distributions (where applicable) and using a 
preferred fragmentation length scale to calculate the core 
mass. We choose to use such methods due to the random 
nature of molecular cloud properties. This allows us to 
statistically determine the shape of the CMF for a wide 
range of randomly chosen cr„ — fiQ pairs. 

2.1. Physical Model 

We consider the formation of cores and the resulting 
CMF within ionized, isothermal, interstellar molecular 
clouds. These clouds are modelled as planar sheets with 
infinite extent in the x- and y- directions and a local ver- 
tical half thickness Z. The nonaxisymmetric equations 
and formulations of our assumed model have been de- 
scribed in detail i n several paper s CCiole k fc Basul[2006l: 
iBasu et al.ll20Q9allbl : iBailev fc Ba su 2012). For this work 
we consider three models: nonmagnetic, flux-frozen mag- 
netic field and a magnetic field with ambipolar diffusion. 



The key ingredient to this analysis is the assumed 
length scale for the core. This length scale for collapse 
can be derived thr ough linear ana. l ysis. The nonaxisym- 
metric equations of lCiolek fc Basul (|2006D and I Basu et al.l 
(J2009a'!b') include the effect of ambipolar diffusion. This 
is quantified by the timescale for collisions between ions 
bound to the magnetic field and free neutral particles. 
This timescale is 

= 1.4 fa^il^) (1) 

\ nii J ni{(Jw)tH2 

Here, is the ion mass, is the number density of 
ions and {(Jw)iH2 is the neutral- ion collision rate. The 
typical atomic and molecular species within a molecular 
cloud are singly ionized Na, Mg and HCO which have 
a mass of 25 amu. Assuming collisions between H2 and 
HC O"*", the neutral- ion collisio n rate is 1.69 x 10~^ cm"^ 
s~^ (jMcDaniel fc Masonlll973l ). Collisions between neu- 
trals and ions transfer information about the magnetic 
field to the neutral particles. The threshold for whether 
a region of a molecular cloud is stable or unstable to col- 
lapse is given by the mass-to-flux ratio of the background 
reference state 

Mo^2.Gi/2g^, (2) 

i>ref 

where (27rG^/^)~^ is the critical mass-to-flux ratio for 
gravitational collapse in the adopted model and 5ref is 
the magnetic field strength of the reference state. Re- 
gions with /xp < 1 are defined as subcritical, regions with 
/xq > 1 are defined to be supercritical and regions with 
/xq ^ 1 are transcritical. 

A dispersion relation for the governing magnetohy- 
drod ynamic equations can be found via linear analysis 
(Ci dek fc Basul [2OO6I : [Basu et al.l [20091 IBailev fc Basul 
2OI2I) . Here we fo llow the analysis as described in 
Bailev fc Basul ()2012|) . For a model with ambipolar dif- 
fusion, the resulting dispersion relation is 

{w + ie){Lu^ ~ C^s^ok^ + 27rGcr„,ofc) 

= uj{2TTGan,oktio^ + k^Vlo) (3) 

where 

6' = T„,o(27rGcr„,ofcM(7^ + fc^V'jo). (4) 

Here, w is the angular frequency of the perturbations, 
Tni,o is the initial neutral-ion collision time, k is the 
wavenumber in the z-direction, Va.o is the Alfven speed, 
where ^ 

F| = = 2^Ga„,oMo '^0, (5) 

47rp„,o 

Zq is the initial half-thickness of the sheet, and Gcff.o is 
the local effective sound speed, such that 

^2 _ 2 [3Pcxt + {7T/2)Gal„] 2 
Coff,o - 2^"-°[P,., + (7r/2)G<o]2^- 

Here, Cg — (fc^T/mn)^/^ is the isothermal sound speed, 
fcs is the Boltzmann constant, T is the temperature in 
Kelvins and m„ is the mean mass of a neutral particle 
(m„ = 2.33 amu). For this analysis, we assume a tem- 
perature T = 10 K and a normalized external pressure 

Pext = 2Pext/7rG(72_o = 0.1. 
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Figure 1. Wavelength with minimum growth time as a function 
of initial mass-to-fiux ratio. Displayed curves are for T„ifi/to = 
(solid curve, flux freezing) and T„ifi/to = 0.2 (dotted curve). 
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Figure 2. Model lognormal column density distribution. 



In the hmit of flux freezing, T„i^o ~^ 0, which gives the 
reduced dispersion relation 

w2 + 2^Ga„,ofc(l - Mo"') - fc2(C^cff,o + Vlo) = 0. (7) 

The gravitationally unstable mode corresponds to one of 
the roots of < and occurs for > 1- The growth 
time for this mode can be written as 



A 

" 2^[Ga„,o(l-A*o')(A-AAfs)]i/2 
for A > Ajv/s, where 



Ams — 



Gct„,o(1 - Mo"') 



(8) 



(9) 



The length scale corresponding to the minimum growth 
time is Ag,m = 2Ams- This is the length scale used to 
produce our synCMFs for models with flux freezing. The 
variation of this length scale as a function of fiQ is shown 
by the solid line in Figure [1] For the case with no mag- 
netic field, Equation |9] reduces down to the thin disk 
equivalent of the Jeans length. 



A,7 = 



G(T„_n 



(10) 



Again, the length scale corresponding to the minimum 
growth time is Xg,m,j = 2Aj, which is the scale used in 
our nonmagnetic model. 



The addition of ambipolar diffusion complicates the 
process somewhat. In these cases, the gravitationally 
unstable mode corresponds to one of the roots of the full 
dispersion relation (Equation [3]). However since it is a 
cubic function, there is no simple expression to describe 
these roots. Therefore, each length scale is computed nu- 
merically. The value of this length scale is related to the 
degree of ambipolar diffusion i.e., the degree of ioniza- 
tion within the cloud, and the mass-to-flux ratio of the 
region. Previous studies show that the ionization frac- 
tion within a mole c ular cloud resembles a step function 
(jRuffle et al.l 119981: iBailev fc Basul I2012D such that the 
outer layers are highly ionized due to UV photoioniza- 
tion while ionization of denser inner regions is primarily 
due to cosmic rays. For this study, we choose to fix 
the neutral-ion collision time to the dimensionless value 
Tni,o/to = "^T^GanfiTnifi/cs = 0.2 ; a value typical of the 
denser inner region s where most cores are likely to form 
(|Basu et al.l l2009al and references within). This corre- 
sponds to an ionization fraction Xi — 5.2 x 10^® at a neu- 
tral column density (t„.o — 0.023 g cm"^. Figure [T] (dot- 
ted line) shows the relation between the collapse length 
scale and the mass-to-flux-ratio for this neutral-ion col- 
lision time. By fixing the neutral-ion collision time, our 
ambipolar diffusion models have only two free parame- 
ters, the column density and mass-to-flux ratio distribu- 
tions. Our choices for these two parameters are discussed 
in the following sections. 
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Figure 3. Model mass-to-fiux distributions for flux freezing mod- 
els. Left: Broad Lognormal Distribution (FF2). Right: Narrow 
Lognormal distribution (FF3). 

2.2. Column Density Distribution 

A survey of column density (T„ distributions within 
various molecular clouds shows that they generally ex- 
hibit log-normal distributions either with or without a 
high density tail (jKainulainen et al]|20090 . Correlation 
of these different shapes with the conditions within the 
clouds suggest that regions with a pure lognormal distri- 
bution tend to be quiescent while those with high density 
tails show signs of active star formation. 

Since the aim of this paper is to investigate the shape 
of the core mass function as an initial condition for 
star formation, we choose a simple lognormal distribu- 
tion as shown in Figure [H This plot shows the distri- 
bution as a function of both the column density ((t„, 
lower axis) and the visual extincti on (A„, u pper axis). 
Following the prescription of Pine da et al.l ()201Cl[ ). the 
conversion from visual extinction to column density is 
achieved by combin ing the ratio of H2 column density to 
color excess ( Bohli n et al. 119781 ) with the total selective 
extinction (,Whittetl 12003 ) to yield a conversion factor 
N{H2) = 9.35 X W°Ay cm-2 mag^^. Although this 
conversion is specifically for H2, the abundance ratio of 
CO to H2 is ~ 10"** and other molecular contributions 
are even smaller, so they do not add significantly to the 
number density of H2 ■ Therefore we assume this number 
density is representative of all species. Assuming a mean 
molecular weight of 2.33 amu, this translates into a mass 
column density conversion of the form 

On = 3.638 X IQ^^A^ g cm'^ mag^^ (11) 

The variance and mean (cr^ and /i) of this distribu- 
tion were chosen based upon observational information. 
Previous studies of molecular clouds show visual ex- 
tinction thresholds for core and star formation to be 
on the order of A„ = 5 mag pohnsto ne et al.! 120041: 
Kirk et "aI]|2006D and A„ = 8 m ag (see P^ohnston e et al.l 
20041 : iFroebrich fc R owTcs '20101 among others) respec- 
tively. As such, we adopted a mean visual extinction 
value of 8 magnitudes for our lognormal density distri- 
bution. The variance reflects the typical width of the 
lognormal fits to cloud density functions presented by 
iKainulainen et al.l (|2009[) . 

2.3. Mass-to-Flux Ratio Distributions 

Although density/ visual extinction maps are fairly 
commonplace, measurements of magnetic field strengths 
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Figure 4. Model mass-to-fiux distributions for ambipolar diffu- 
sion models. Left: Broad Lognormal Distribution (AD4), Right: 
Narrow Lognormal (ADS). 



within molecular clouds are difficult to obtain. Due 
to limitations in techniques and resolution, studies 
of magnetic fields wi thin c louds are generally on a 
more global s cale (se e [Crutchcr 1999: Heilcs & Tr olan dl 
'2004; Troland fe Crutchcr 2008; Falgaro ne et al.l 120081 
iCrutcher et al.i i2010. : ,Chapman et al.. ,201ll among oth- 
ers) which does not give much insight into the exact na- 
ture of /zo within denser small scale regions. Therefore, 
the mass-to-flux ratio of speciflc regions are not generally 
known, let alone a distribution over an entire cloud. Re- 
cent simulations of cloud forma tion with magnetic flelds 
()Vazauez-Semadeni et al.ll2011[ ) show that the mass-to- 
flux ratio distribution seems to exhibit a lognormal 
shape. On the other hand, analysis o f the likelihood 
of diff erent magnetic field distributions (jCrutcher et al.l 
120101 ) show that the magnetic field strengths for various 
regions (HI diffuse clouds, OH dark clouds, etc) exhibit a 
uniform distribution ranging from very small values up to 
a maximu m value. This seems to disagree with the sim- 
ulations of iVazguez-Semadeni et all ()2011[ ). With these 
results in mind, we choose to explore both options (i.e., 
uniform and lognormal distributions). 

As shown by the l inear analysis results presented in 
iBailev fc Basui (|2012[ ) and Figure [TJ the length scale for 
collapse is dependent on the value of the mass-to-fiux 
ratio. The value of /xg is selected from a predefined 
distribution that is independent of the distribution of 
iT„. This implies that the magnetic field strength is not 
constant and varies according to the choices of cr„ and 
^iQ. The independent sampling of values of an and Ho 
does not then allow for any systematic dependence of 
one quantity on the other. We believe this is an accept- 
able first approximation since the initial conditions of 
the mass-to-fiux ratio distribution in a molecular cloud 
are poorly constrained. We test several possible /xq dis- 
tributions in an attempt to determine if the shape of 
an observed CMF could reveal information about the 
underlying mass-to-flux ratio distribution. We consider 
both uniform and lognormal distributions. Figures [3] & |4] 
show the adopted lognormal mass-to-flux ratio distribu- 
tions for the flux freezing and ambipolar diffusion mod- 
els respectively. Specifically, all distributions sample the 
transcritical peak in fragmentation scale, Xg^m (see Fig- 
ure [T]). The properties of all fiQ distributions considered 
are given in Table [T] 
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Table 1 

Model Parameters 



Model Name 


fiQ Distribution 


Mean (fi) 


Variance (cr'^) 


fiQ Range 


Flux Frozen Models 


FFl 


Uniform 






1.0 - 3.0 


FF2 


Broad lognormal 


0.01 


1.0 


1.0 - 10 


FF3 


Narrow lognormal 


0.01 


0.01 


1.0 - 1.5 


Ambipolar Diffusion Models 


ADl 


Subcritical Uniform 






0.1 - 1.0 


AD2 


Supercritical Uniform 






1.0 - 3.0 


ADS 


Uniform 






0.7 - 3.0 


AD4 


Broad lognormal 


0.01 


1.0 


0.3 - 10 


ADS 


Narrow lognormal 


0.01 


0.01 


0.6 - 1.5 




-2-10123 .2-10123 
Log(Mass) (M^J LogCMass) (M^ 

Figure 5. Synthetic core mass function for a non magnetic cloud. 
Left: Total core mass function. Right: Contributions to the core 
mass function from cores with Ay < 8 mag (dashed line) and cores 
with Ay > 8 mag (dotted line). 



2.4. Producing Synthetic Core Mass Functions 

To produce a synthetic CMF, we randomly sample the 
column density distribution for the nonmagnetic case and 
both the column density and mass-to-flux ratio distribu- 
tions for the magnetic cases. These values are then used 
to flnd the preferred length scale for collapse from the 
linear analysis. Finally, the mass is determined by mul- 
tiplying the column density by the square of the cor- 
responding length scale. By randomly sampling each 
model distribution 10^ times, a synthetic CMF is pro- 
duced. 

3. MODELS AND RESULTS 

Our analysis covers several different mass-to-flux ratio 
distributions and assumptions about the neutral-ion col- 
lision time and column density distribution. As stated 
earlier, the column density distribution is the same for all 
models (see Figure [5]) and the neutral-ion collision time 
for the ambipolar diffusion models is set to a normalized 
value, Tni o/to = 0.2. In addition to the models listed 
in Table [Tl we also present a nonmagnetic (NM) fidu- 
cial case. The following subsections present the results 
for each model individually. An in depth comparison be- 
tween all the models and implications regarding observed 
CMFs will be discussed in Sections 13.41 & |3] respectively. 

3.1. Non-Magnetic Model 

The nonmagnetic model serves as a baseline for our in- 
vestigation. The left panel of Figure[5]shows the resulting 
core mass function from this technique. As discussed in 
Section 12.21 we choose the peak of our density distribu- 
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Figure 6. Synthetic core mass functions for a flux frozen magnetic 
cloud assuming a uniformly distributed mass-to-flux ratio (FFl). 
Left: Total core mass function. Right: Contributions to the core 
mass function from cores with Ay < 8 mag (dashed line) and cores 
with Ay > 8 mag (dotted line). 

tion to correspond to the apparent visual extinction 
threshold for the creation of star forming cores; ^ 8 
magnitudes. The right hand panel of Figure [5] shows the 
contributions from high density gas (A^ > 8 mag, dot- 
ted line) and low density gas {Ay < 8 mag, dashed line). 
As expected from the Jeans theory, the core mass dis- 
tribution mimics the column density distribution, with 
high mass cores formed from low density gas and low 
mass cores formed from high density gas. The distri- 
bution of masses for this model peaks at a value of 
log(M/MQ) = 0.4 or -M" 2.5Mp, which is consis tent 
with observations (jNutter fc Ward-Thompsonll2007f) . 

3.2. Flux Frozen Magnetic Model 

A main aim of this paper is to show the effect of a mag- 
netic field on the CMF. A flux frozen field represents the 
simplest case. Such a scenario arises in highly ionized re- 
gions where frequent collisions between ions and neutral 
particles would ensure perfect coupling to the magnetic 
field. Figures |6]|8] show the resulting synthetic core mass 
function for the three models FFl, FF2, and FF3 respec- 
tively. Under the assumption of a uniform mass-to-flux 
ratio distribution (FFl), the resultant CMF (Figure |6l 
left) exhibits a narrow peak with a distinct high mass 
tail. The right hand panel of Figure [6] again shows the 
contributions to the CMF from the two column density 
regimes {Ay < 8 mag (dashed line) and Ay > 8 mag (dot- 
ted line)). This composite plot shows that like the NM 
case, and in line with the Jeans theory, the low density 
gas forms high mass cores, while high density gas forms 
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Figure 7. Synthetic core mass functions for a fiux frozen mag- 
netic cloud assuming a broad, lognormal mass- to- fiux ratio (FF2). 
Panels depict the same curves as Figure [G] 
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Figure 8. Synthetic core mass functions for a fiux frozen mag- 
netic cloud assuming a narrow, lognormal mass-to- flux ratio (FF3). 
Panels depict the same curves as Figure [S] 
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Figure 9. Synthetic core mass functions for a magnetic cloud 
including the effects of ambipolar diffusion assuming a uniform 
subcritical distributed mass-to-flux ratio (ADl). Panels depict the 
same curves as Figure [G] 




Log(Mass) (M^J Log(Mass) (M^J 

Figure 10. Synthetic core mass functions for a magnetic cloud 
including the effects of ambipolar diffusion assuming a uniform 
supercritical distributed mass-to-flux ratio (AD2). Panels depict 
the same curves as Figure |6] 



low mass cores. However, unlike the Jeans theory and 
nonmagnetic case, we see that with the addition of a 
magnetic field, the high density gas also contributes to the 
formation of high mass cores, albeit to a lesser extent. 
Compared to the NM case, the peak of this core mass 
function is shifted to M ~ 10° '^ Mq ~ 5.0 Mq. On the 
right hand side of this peak, the trend can be described 
by two distinct slopes. For 0.7 < log(M/M0) < 1.2, 
a — 0.8 while for log(M/MQ) > 1.2 the slope becomes 
shallower; a ~ 0.6. Neither of these values correspond 
to the typical Salpeter and observational values. This 
discrepancy will be discussed further in Sectional 

The formation of the high mass tail is due to the re- 
lationship between /xq and A as defined by Equation [9l 
For /zo — cr„ pairs which have mass-to-fiux ratios closer to 
the critical value (/zq = 1), the corresponding length is 
up to 23 times larger than the thermal Jeans length for 
the same column density (see Figure [!}. This increase in 
length scale has a direct effect on the mass of the core 
that is formed. Conversely, the low mass distribution is 
formed by fiQ — f7„ pairs that have low column density 
and mass-to-flux values that are closer to the other limit 
(/io = 3), where A is only about 1.5 times larger than the 
thermal length scale. 

Figure [7] shows the resulting synCMF for a broad log- 
normal /-to distribution (FF2). The two panels again 
show the total and composite CMFs as described above. 



This distribution results in a CMF that is similar to that 
of model FFl (Figure (5]), with a few minor differences. 
First, the high mass tail exhibits a steeper slope that 
results in a more pronounced peak region. Second, the 
peak of the mass function has shifted to a slightly smaller 
value of M ~ 10°-^ Mq = 3.16 Mg. As before, the trend 
of the high mass side can be described by two distinct 
slopes. For 0.5 < log(M/M0) < 1.0, a = 1.31 while for 
log(-M/M0) > 1.0 the slope becomes shallower; a — 0.63. 

Figure [5] shows the resulting synCMF for a narrow log- 
normal fiQ distribution (FF3). Unlike the previous two 
models, this one does not exhibit a narrow log-normal 
type peak, but rather shows a broad peak that leads 
directly into a high mass tail. As a result, the post peak 
trend for this model can be described by a single slope, 
a — 0.44. Also, note that the function itself has been 
shifted toward higher masses as compared to the other 
two flux frozen models. As such, this CMF peaks at 
M ~ lO^-^ M0 - 20 Mq. This shift in the mass range is 
due entirely to the narrow peak distribution of the mass- 
to-flux ratio; all of the chosen mass-to-flux ratios result 
in length scales that are ~ 6 — 23 times larger than the 
thermal length scale (see Figure[T]) and therefore, the low 
mass cores that are formed in the other two models are 
absent in this model. Overall, as shown by all three 
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models, the effect of adding a flux-frozen field is the ap- 
pearance of a broad shallow tail at the high mass end of 
the core mass function. 



Log(Mass) (M^J Log(Mass) (M^J Log{Mas5) {M,J 

Figure 11. Synthetic core mass functions for a magnetic cloud including the effects of ambipolar diffusion assuming a uniformly distributed 
mass-to-flux ratio (ADS). Left: Total core mass function. Middle: Contributions to the core mass function from cores with < 8 mag 
(dashed line) and cores with Ay > 8 mag (dotted line). Right: Contributions to the core mass function from cores with fio < I (dashed 
line) and cores with > 1 (dotted line). 

above. Specifically, this CMF shows the same peaked na- 
ture with high mass tail as the flux frozen model, however 
this tail abruptly declines at about 100 M©. This trun- 
cation makes the over all shape of the CMF resemble the 
nonmagnetic case, albeit broader, with the beginnings of 
a "shoulder" feature between 10 and 100 Mq. Looking 
at the composite column density CMF (Figure [101 right), 
we see that the lowest and highest mass cores are formed 
by the highest and lowest density gas respectively, while 
the middle has contributions from both density regimes. 
The peak of the mass function for this model occurs at 
about log(M/M0) = 0.7. 

Model AD3 assumes a uniform mass-to-fiux ratio dis- 
tribution that samples the peak of the A versus fiQ graph 
(see Figure [T]). The resuhing CMF (Figure [TTl left) is 
very similar to the one produced by AD2. Looking at 
the contributions from the low and high column density 
gas, we again see that the lowest and highest mass cores 
are formed by the highest and lowest density gas respec- 
tively while the middle range has contributions from both 
density regimes. 
The right panel of Figure [TT] shows the contributions 



3.3. Ambipolar Diffusion Magnetic Model 

In the previous section we looked at the effect of a 
simple flux-frozen field on the shapes of the resulting 
CMF(s). Here we look at how the addition of neutral- 
ion slip via ambipolar diffusion affects the shape of the 
CMF. As discussed above, we have fixed the normalized 
neutral-ion collision time to Tni,o/to — 0.2. This implies 
a high degree of ambipolar diffusion and therefore less 
frequent collisions between the neutrals and ions. Such a 
situation would occur in the inner regions of a molecular 
cloud where the main ionization mechanism is cosmic 
rays. 

Figures l9lfT3l show the resulting synCMFs for all five 
mass-to-fiux ratio distributions respectively. To estab- 
lish how the sub- and supercritical regions of the mass- 
to-flux ratio affect the shape of the CMF, we start our 
analysis by presenting two cases that isolate each regime. 
Figures ^ & [TU] show the resulting synCMFs for the sub- 
critical and supercritical uniform mass-to-fiux ratio dis- 
tributions (ADl and AD2) respectively. The two panels 
show the total and constituent core mass functions as 
described in the previous section. 

Focusing on model ADl, Figure [HI the left panel shows 
that the core mass function is very similar to the non- 
magnetic model (see Figure [H left). This is due to the 
fact that the curve on the subcritical side of Figure [1] con- 
verges to the nonmagnetic limit faster than in the trans- 
and supercritical regions. Upon closer comparison, ADl 
peaks at the approximately the same value as NM, how- 
ever the density composite CMF (Figure [HI right) reveals 
differences between these two models. Unlike the non- 
magnetic model, ADl shows evidence that a portion of 
the high column density gas goes toward forming high 
mass cores (Figure [9l right). 

Figure [TUj shows the resulting synCMF under the as- 
sumption of a uniform supercritical distribution (AD2). 
The left panel shows that the total CMF is a hybrid be- 
tween the nonmagnetic and flux-frozen models presented 



from the subcritical (/io < 1, dashed line) and super- 
critical (/io > 1, dotted line) gas. We see that the total 
synCMF for ADS (Figure [11] left) is a combination of 
models ADl and AD2. Specifically, we see that the ma- 
jority of the cores are formed from supercritical gas, while 
the subcritical gas yields a minor contribution to the pop- 
ulation of low mass cores. By mentally combining the 
middle and right hand plots in Figure [11] one can deter- 
mine that the highest mass cores are formed by supercrit- 
ical gas and fall into the non-star-forming regime while 
low-mass cores are formed by both supercritical and sub- 
critical gas, and fall into both the star-forming and non- 
star-forming regimes. The peak of the mass function for 
this model occurs at about log(M/M0) = 0.7 and the 
average slope of the high mass 'tail' is a = 1.42. 

Finally, Figures [12] & [El show the resulting synCMFs 
for the two lognormal /io distributions, AD4 and ADS, 
respectively. The broad lognormal distribution (AD4) is 
similar to models AD2 and ADS, however this model 
shows a more distinct 'peak' and 'shoulder' region as 
compared to the other two. Looking at the composite 
mass-to-fiux ratio plot (Figure [121 right) we see that 
the peak region is mainly formed by subcritical gas 
while the shoulder region is formed mainly by contri- 
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Figure 12. Synthetic core mass functions for a magnetic cloud including the effe cts of ambipolar diffusion assuming a broad, lognormal 
mass-to-flux ratio distribution (AD4). Panels depict the same curves as Figure [TT] 
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Figure 13. Synthetic core mass functions for a magnetic cloud including the effec ts of ambipolar diffusion assuming a narrow, lognormal 
mass-to-flux ratio distribution (ADS). Panels depict the same curves as Figure [TT] 



butions from supercritical gas. This model peaks at 
M = 10°-'^ M/Mq ~ 5.0 M©, and the average slope 
of the high mass tail is a = 1.18. Switching to the nar- 
row lognormal distribution (Figure I13p , we see that this 
model results in a double peaked function. Examination 
of the composite plots show that the low mass peak is 
formed by the subcritical material while the second peak 
is formed by supercritical material. These peaks occur 
at log(M/M0) - 0.7 and log(M/M0) - 1.5 respectively. 
The formation of the high mass peak is due to the ex- 
tremely narrow mass-to-flux ratio distribution. It picks 
out only large length scales from the peak of the X — Ho 
curve (with Tnifl/to = 0.2) in Figure [II 



3.4. Assessment of Synthetic Core Mass Functions 

The previous subsections presented the overall results 
and features of each of the models. Within these results 
we found three main features that changed between the 
different models. These are the overall shape of the core 
mass function, the location of the peak(s) and the slope 
of the high mass tail (if it exists) . Here we discuss these 
three features across all models. 



3.4.1. Shape 

Within the nine models presented, there were three dis- 
tinct recurring shapes; pure lognormal as represented by 
the NM and ADl models, lognormal peak with a shoul- 
der as represented by AD2, ADS, AD4 and AD5, and 



the lognormal peak with high mass tail as represented 
by FFl, FF2, and FF3. The appearance of these shapes 
are directly connected to the state of the magnetic field 
in the region. In the absence of a magnetic field, the 
CMF is a pure lognormal function. This shape is also 
observed in model ADl. As mentioned earlier, the rea- 
son that this AD model shows such a shape while the 
other ones do not is due to the shape of the A — /xg curve 
on the subcritical side of Figure [T] the curve asymptotes 
to the nonmagnetic limit faster on that side than on the 
supercritical side. Therefore one would expect a model 
with only subcritical mass-to-flux values to look similar 
to the nonmagnetic model, but with a slight broaden- 
ing due to a narrow region of mass-to-flux ratios with A 
larger than the non-magnetic limit. 

For models with an increasing supercritical regime, the 
broadening becomes more pronounced as a shoulder de- 
velops. This shoulder is due to an increase in higher 
mass cores that are the product of the larger length scales 
picked out by the supercritical mass-to-flux ratios. The 
extent of the shoulder depends on the mass-to-flux ra- 
tio distributions. For uniform distributions, the CMF is 
narrower with a less defined shoulder region, while for 
a broad lognormal distribution, the shoulder region is 
much broader and distinct. Finally, the appearance of 
the double peaked CMF in AD5 is an example of an ex- 
treme shoulder. This second peak is due solely to the 
extremely narrow mass-to-flux ratio range used in this 
model. This preferentially picks out only mass-to-fiux 
ratios with length scales much larger than the nonmag- 
netic model. 

The appearance of the pure high mass tail is entirely 
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Figure 14. Small sample core mass functions for three models: NM (top row), FFl (middle row), and ADl (bottom row). Number of 
cores for each panel indicated in the top left hand corner of each plot. 
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Figure 15. Bin size comparison for small sample core mass functions. Panels show the effect of the bin size on the resulting curve for 
Alog(Af/Af0) = 0.25 bins (top row) and Alog(M/M0) = 0.1 bins (bottom row). Model used in all panels is ADS. 

a product of flux-freezing. This is due to the asymptotic 
nature of the flux-frozen curve as it nears the critical 
mass-to- flux ratio (see Figure [T]). This allows for trans- 
critical mass-to-flux ratio values to produce much larger 
masses for the same column density. 

3.4.2. Peak Location 

The location of the CMF peak depends on the distribu- 
tion of the mass-to-flux ratio. The location of the peak in 
the nonmagnetic case, which occurs at log(Af/M0) = 0.4 
(M ~ 2.5 M0) serves as the comparison point. For 



magnetic models, the location of the peak was gener- 
ally larger than this value as long as the mass-to-flux 
ratio distribution was uniform with some contribution 
from the supercritical regime (see models FFl, AD2, and 
ADS). Model ADl, although also assuming a uniform 
mass-to-flux ratio distribution, exhibits a similar peak 
value to NM due to the exclusion of supercritical mass- 
to-flux ratio values. When considering the lognormal 
mass-to-flux ratio distributions, we find that the peak 
location is dependent on the width of the distribution. 
Specifically, broader distributions exhibit values closer 
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to the NM peak value, while narrower distributions ex- 
hibit peak values that are higher than the nonmagnetic 
case. Model AD5 is an anomaly and does not fit within 
these trends given that it exhibits two peaks. 

3.4.3. High Mass Slope 

As alluded to earlier, the shape and extent of the high 
mass slope was found to be variable and connected to 
the influence of the magnetic field. Specifically, the ap- 
pearance of the 'shoulder' feature is directly connected 
to the presence of ambipolar diffusion. The degree of the 
shoulder in the ambipolar diffusion models was found to 
be dependent on the range of allowed mass-to-flux ra- 
tio values. Overall, these differences in shapes result in 
a wide range of slopes. For the flux-frozen models, the 
slopes were as steep as a = 1.31 in the case of FF2, 
and as shallow as a = 0.44 in the case of FF3. For the 
ambipolar diffusion models, the average high mass slope 
ranges between a = 1.18 and a — 1.42. Although some 
of these slopes arc consistent with the Salpeter value, 
a = 1.35 (Salpeter 1955), others are significantly differ- 
ent. Further analysis of this discrepancy is given in the 
following section. 

4. SCALING TO OBSERVATIONS 

Unlike our synCMFs, typical observational CMFs usu- 
ally contain on the order of 200 cores, not 10^. Therefore, 
to make our analysis relevant for typical observed CMFs, 
we must scale back our sample sizes to those typically ob- 
served. The following two sections explore the effect of 
two observational constraints, sample size and bin size, 
on the shape and slope of observed CMFs. 

4.1. Effect of Sample Size 

To test the effect of the sample size on the resultant 
CMF, we scaled three synCMFs (NM, FFl, and AD3) 
down to plausible observational sample sizes (100, 200, 
300, 400, and 500 cores). Figure [U] shows the result- 
ing synCMFs for each of the fifteen cases. In addition 
to scaling the sample size, we have also truncated the 
mass range considered to one more typically found in ob- 
served CMFs (-1.0 < log(M/Mo) < 1.3). Under these 
scaled conditions, we see that the nonmagnetic CMFs 
still maintain the overall shape exhibited by the full sam- 
ple curve (Figure [5]) , however the two magnetic cases are 
fairly different. The high mass tail and truncated shoul- 
der features present in the full sample curves for FFl 
and AD3 respectively are no longer quite as distinct at 
these sample sizes. For a definitive difference between 
the ambipolar diffusion and fiux-frozen cases, observa- 
tions would have to extend up to objects with masses 
between 10^ and 10'^ solar masses. Therefore, on typical 
observational scales, conclusions about the nature of the 
magnetic field from the shape of the CMF are possible, 
but highly uncertain. 

4.2. Effect of Bin Size 

Constructing histograms for the purposes of deter- 
mining a CMF requires binning data into predeter- 
mined mass bins. For the above synCMFs, we used 
Alog(M/MQ) = 0.1 size bins. Variations in the bin 
size acts to change the resolution of the resulting curve; 
smaller bins yield more detail while larger bins show only 
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Figure 16. Average slope as a function of sample size. Sym- 
bols represent the derived slopes for the two models and two bin 
sizes: ADS, Alog(M/M0) = 0.1 (squares), ADS, Alog(M/M0) = 
0.25 (circles), FFl, Alog(M/M©) = 0.1 (diamonds), and FFl, 
Alog(M/MQ) = 0.25 (triangles). Average slopes computed over 
a minimum of 2000 samples. Open symbols indicate the slopes of 
the full sample size. 

the broad strokes. To determine the effect of the bin size 
on the resulting CMF, we re-binned the histograms for 
AD3 in Figure [Til (bottom row) using Alog(M/M0) = 
0.25 bins. Figure [TS] shows the comparison of the original 
bin size (Alog(M/MQ) = 0.1) to the new bin size (top 
row). As expected, with the larger bin size, the detail 
becomes smeared out, resulting in an average curve. 



4.3. Effect of sample size and bin size on CMF slopes 

The main piece of data generally extracted from a 
CMF is the slope of the high mass tail. This information 
is then used to compare different regions to each other, 
and to the initial mass function (IMF) in an attempt 
to determine the true nature of star formation and the 
possible relation between the CMF and IMF. However, 
as discussed above, the sample size and bin size have a 
profound effect on the shape of the curve. This effect 
translates over to the derived slopes. To determine the 
extent of this effect, we generate multiple small sample 
CMFs (2000-I-) for each sample size and compute the 
average slope. Figure [12] shows the results of this anal- 
ysis for models FFl and AD3 for both mass bin sizes. 
The filled symbols show the average slope for each of the 
models while the open symbols depict the slope of the 
full sample (10^). Tests with larger numbers of samples 
for each sample size showed differences in the average 
slope of up to 0.01, which is encompassed in the size of 
the symbols. 

As shown in Figure[T6l the size of the bin clearly affects 
the average slope. The larger bin size yields slopes that 
are steeper than the Salpeter slope, while the smaller 
bin size shows an overall shallower average slope. The 
size of the sample also effects the slope. Smaller samples 
generally result in steeper slopes than those derived using 
the full sample. Furthering this analysis we look at both 
the minimum and maximum slopes calculated for 
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Figure 17. Histogram of computed slopes for different models 
and sample sizes assuming Alog(M/M0) = 0.1 bins. Top row: 
ADS with 100 cores per sample (left) and 500 cores per sample 
(right). Bottom row: Same as the top but for FFl. 




Figure 18. Histogram of computed slopes assuming a CMF with 
100 cores and Alog(M/MQ) = 0.25 bins. Left: ADS, Right: FFl. 



each filled point in Figure (THl as well as the distribution 
of slopes. Figure [17] shows the distribution of slopes for 
four of the points on Figure [16] as indicated (Top row: 
ADS, Bottom Row, FFl. Left column: 100 cores. Right 
column: 500 cores) assuming a CMF constructed with 
Alog(M/M0) = 0.1 bins. Ah four cases show that the 
preferred slope value is close to the average slope value. 
The maximum and minimum computed slopes exhibit a 
very wide range for the small sample sizes (i.e., a = —0.1 
to a = 11 for ADS, 0.1, 100 cores) while the larger sample 
sizes exhibit a smaller maximum- minimum range (i.e., 
a = 0.S4 to a = 1.79 for AD3, 0.1, 500 cores). This 
decrease in the slope variance is evident when comparing 
the left column to those in the right column in Figure [17] 
From these plots, we conclude that although there can 
be a wide variance in possible slope values, the preferred 
slope value is in general smaller than the typical Salpeter 
value, a — 1.35, and the range of slopes decreases as the 
number of samples increases. Figure [T5I shows the distri- 
bution of slopes for two CMFs (left: ADS, right: FFl) as- 
suming a 100 core sample size and Alog(Af/Af0) — 0.25 
bins. Comparing to the left hand plots in Figure 1171 
see that the larger bin size results in a bimodal distribu- 
tion with the peaks occurring at a « 2.0 and a « 1.5 
for ADS and a « 1.9 and a « 1.45 for FFl. The resuh 
of this bimodal distribution is to shift the average slope 
values to smaller values than the dominant peak. This 



is particularly evident in Figure [TC] in the trend of slopes 
for the smallest sample sizes of the blue triangles (FF, 
0.25). Further analysis of the effect of the original col- 
umn density distribution on the variance and mean of the 
resulting slope histogram showed that a larger variance 
in the column density distribution shifts the mean in the 
slopes to smaller values (< 1) while a smaller variance 
results in a larger mean value, a ^ 1.S5. 

5. DISCUSSION 

Our analysis shows that the shape of the CMF is highly 
dependent on the magnetic field strength and neutral- 
ion coupling within the cloud. Specifically, a fiux-frozen 
magnetic field broadens the nonmagnetic lognormal dis- 
tribution to have a significant power-law high mass tail, 
though it is much shallower than the Salpeter value. 
When ambipolar diffusion is taken into account, there 
is an intermediate mass tail and a high mass cutoff. The 
extent of all these features are dependent on the range 
o f mass-to-flux ratio val u es in the initial cloud. 

iKunz fc Mouschoviad (|2009l ) carried out a more fo- 
cused study of the effect of magnetic fields and ambipolar 
diffusion in creating a broad CMF. Their model explored 
only the subcritical portion of the fragmentation scales 
seen in Figure [1] Furthermore, they assumed a uniform 
distribution of subcritical mass-to-flux ratios and effec- 
tively a fixed Jeans mass in order to generate their mass 
distribution. 

The low-mass tail in their distribution originates in the 
assumption that the subcritical clouds ultimately form 
dense cores with masses that are scaled by /ig for subcrit- 
ic al values of un ■ This i s beca use numerical simulations 
of iBasu Sz MouschoviasI (|1995l ) show that only an inner 
region where the mass-to-flux ratio exceeds the critical 
value undergoes rapid collapse. We do not make that as- 
sumption in this study, since cores that form by ambipo- 
lar drift have an appearance that is similar to those that 
are forming by a more rapid gravitationally-dominated 
process (see Basu et al. 2009b). Since the resultant CMF 
in our model is generated from an underlying lognormal 
function, it has an intrinsic peak even when binned in lin- 
ear mass bins. An advantage of the KM09 model is that 
they do not need to assume an underlying lognormal dis- 
tribution to obtain a lognormal-like CMF, however that 
CMF is peaked only when binned in log mass. 

Upon scaling our models down to observational sam- 
ple sizes and ranges, we found that the distinction be- 
tween the different models is lost within typical obser- 
vational mass ranges and therefore no information re- 
garding the magnetic field can be reliably gleaned from 
the shape of the observed CMFs. Further to this, anal- 
ysis of the slopes for each of the sample sizes showed 
that the smaller samples sizes result in slopes that are 
1.1 — 1.4 times larger than the slope derived from the 
full sample, while the derived slopes for the larger bin- 
size are ~ l.S — 2.0 times larger than the corresponding 
smaller binsize slope measurements. Although we have 
taken care to scale our analysis down to those typically 
used in observations, the question still remains as to how 
well our results and conclusions correspond to actual ob- 
servations. A recent study of the CMF for five sepa- 
rate star forming regions (Ophiuch us, Taurus, Perseus , 
Serpens and Orion) performed by iSadavov et al.l (|2010( ) 
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provides the perfect platform for comparison. Looking 
at the core mass distributions for these regions, as ex- 
pected, it is hard to definitively discern any characteris- 
tic features that are indicative of a particular magnetic 
field model. With limited data, it is plausible that the 
CMFs for Ophiuchus, Taurus and Perseus could exhibit 
the indicative shoulder of the ambipolar diffusion mod- 
els, while the full Orion CMF could show evidence of a 
flux frozen fie ld. Looking at t he slo pes of the CMFs for 
these regions, iSadavov et al.l (|201Clf ) showed each region 
gave slope values that are close to the a = 1.35 Salpeter 
slope, within their adopted errors. Comparing their slope 
values to those in Figure most of them would fall 
somewhere in the lower half of the graph in amongst the 
diamonds and squares while the Orion with OMC slope 
would fall in amongst the triangles and circles. How- 
ever looking at the binsize of the observations, all of 
the slopes should be within the triangle/circle regime of 
the graph. Comparing these values to the correspond- 
ing slope histograms (see Figure [TS]) we see that these 
values all fall within the regime of possible slopes. On 
the surface, this seems to be a huge discrepancy between 
our results and observations, however each of these five 
observational slopes represents a single slope within our 
2000-1- values used to derive an average slope. However, 
looking at the range of slopes derived from our analysis, 
these observed slopes fall within this range. As shown 
in Figure fTTl the only way to produce a narrower range 
of slope values is to increase the sample size, which is 
not always possible observationally since the number of 
objects detected depends entirely upon the number of 
objects actually present and the sensitivity of the instru- 
ment. 

Based on our analysis and the above comparison to 
the work by Sadavoy ct al. (2010), we argue that the ob- 
served CMFs are extremely statistically limited, both in 
the size of the sample and the number of samples over 
which the slope of the CMF is averaged. Through our 
analysis, we have shown that with larger number statis- 
tics, not only is the measured slope of the CMF much 
different than the typical Salpeter value a = 1.35, but 
highly dependent on the size of the mass bin. In ad- 
dition, the range of individual slope values within the 
set size decreases as the number of cores in the sample 
siz e increases. This is analogous to the results found 
by lElmegreenl (|1999( ) , where although it was determined 
that the most probable value for the IMF slope is the 
Salpeter value, a — 1.35, it is a highly reduced average of 
all possible outcomes. Subseque ntly, we argue tha t based 
on our analysis and the results of lElmegreenI (|1999[ ). there 
seems to be no clear cut correlation between the slope of 
the CMF and the IMF and that the shape and slope of 
the CMF are entirely controlled by the conditions within 
the cloud itself. Since it is unfeasible to claim that all 
clouds exhibit identical conditions, it is therefore unreal- 
istic to expect a universal shape and slope value for all 
star forming regions. 

6. SUMMARY 

We have studied the effect of magnetic fields on the 
formation and properties of the core mass function us- 
ing a combination of the results from linear analysis and 
Monte Carlo methods. In addition, we have studied the 
effects of low number statistics on the slope of the high 



mass tail. Here we summarize the main results of our 
analysis. 

• The synthetic CMFs show that the presence of a 
magnetic field has several effects on the shape of the 
CMF. In general, a magnetic field acts to broaden 
the core mass function compared to the nonmag- 
netic CMF. In addition, the magnetic CMFs ex- 
hibit a high mass tail. The form of this tail de- 
pends on whether the field is flux frozen or allows 
for neutral-ion drift across the field lines. In the for- 
mer case, the tail exhibits a continuous power law 
while in the latter case, the high mass tail truncates 
to form a "shoulder" . 

• The nonmagnetic model shows that the high mass 
cores are formed from low density gas and vice 
versa. Analysis of the contributions of low and high 
density gas to the low and high mass regions of the 
CMF shows that the addition of magnetic fields re- 
sults in additional contributions of high mass cores 
formed from high density gas. 

• Scaling of the synthetic CMFs down to typical ob- 
servational sample sizes and bin sizes show that the 
ability to distinguish between the different models 
is no longer possible for the smallest sample sizes 
(100 cores) and typical bin sizes (Alog(M/M0) = 
0.25). This shows that the current observations of 
core mass functions are statistically limited. 

• Statistical analysis of the derived slope from a large 
sample of synthetic CMFs show that the slope 
of the high mass tail is systematically steeper for 
smaller core sample sizes than for larger sample 
sizes. In addition, the average slope is also system- 
atically steeper for larger bin sizes. 

• Analysis of the minimum, maximum and distribu- 
tion of calculated slopes shows that the most prob- 
able slope does not necessarily correspond to the 
canonical Salpeter value. In addition, the most 
probable slope value becomes shallower as the sam- 
ple size increases. 
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